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We regularize QCD using the combination of higher covariant derivatives and Pauli- 
Villars determinants proposed by Slavnov. It is known that for pure Yang-Mills 
theory the Pauli-Villars determinants generate unphysical logarithmic radiative cor- 
rections at one loop that modify the beta function. Here we prove that when the 
gauge fields are coupled to fermions so that one has QCD, these unphysical correc- 
tions translate into a violation of unitarity. We provide an understanding of this 
by seeing that Slavnov's choice for the Pauli-Villars determinants introduces extra 
propagating degrees of freedom that are responsible for the unitarity breaking. This 
shows that Slavnov's regularization violates unitarity, hence that it should be rejected. 



1 Introduction. The advantages of using gauge invariant regularization methods are well 
known to the quantum field theory practitioner. The problem is that there are not so many 
gauge invariant regularization methods available. The two most popular ones proposed to 
date are probably dimensional regularization |J and the method of higher covariant deriva- 
tives @. Dimensional regularization works well for vector gauge theories, for which the 
algebraic structure is not altered by a change in the number of dimensions of spacetime. 
Unfortunately, when it comes to chiral gauge theories, it is not clear [[| |4j] whether it is pos- 
sible to consistently define dimensional regularization, the reason being that the properties 
of chiral objects depend on the dimensionality of spacetime and this conflicts somehow with 
the ideas behind dimensional regularization. 

The situation for the method of higher covariant derivatives is more confusing. To de- 
scribe it, we will restrict ourselves to QCD, the theory we will be discussing here. As is well 
known, higher covariant derivatives only provide a partial regularization, since they leave 
one-loop divergences unregularized. To achieve full regularization, a second regulator taking 
care of the unregularized one-loop divergences must be introduced. Choosing such a regula- 
tor is not a simple issue, since one would like to pick one that preserves gauge invariance and 
that at the same time does not jeopardize what has been gained at two and higher loops with 
the introduction of higher covariant derivatives. Slavnov || proposed in the seventies to use 
as second regulator a certain combination of gauge invariant Pauli-Villars determinants [see 
eq. (H) for their expression]. Adopting his proposal, one ends up with a hybrid regularization 
that combines higher covariant derivatives with Pauli-Villars determinants of a certain form 
and that keeps the dimension of spacetime at its physical value. We will call this regulariza- 
tion prescription Slavnov's regularization and denote it by SR. Keeping the dimensionality 
of spacetime unchanged and preserving gauge invariance, SR looks like a good starting point 
to formulate a suitable regularization method for chiral gauge theories. Unfortunately, there 
is no agreement as for whether or not SR is a consistent regularization method. On the one 
hand, there are claims that (i) the Pauli-Villars determinants spoil regularization at two and 
higher loops and (ii) that, even at one loop, renormalization is inconsistent with gauge 
invariance [0] . One the other, there are calls [§] rebating these claims. 

To settle down the dispute, and motivated by its potentiality for chiral gauge theories, 
SR was used in ref. || to explicitly regularize and renormalize Yang-Mills theory at one 
loop. There it was proved that the Pauli-Villars determinants on which SR is based generate 
unphysical logarithmic radiative corrections that modify the beta function of the theory at 
one loop, giving for the latter an unphysical value. The purpose of this paper is to show that 
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these unphysical corrections produce a violation of unitarity when the Yang-Mills fields are 
coupled to fermionic matter, so that one has full QCD. 



2 Slavnov's regularization. Let us start by briefly recalling the basics of SR. We are interested 
in iV c -coloured, iVf-flavoured QCD in four dimensions. Not to fall short of rigour in the 
computation of Feynman integrals, we will work in Euclidean space and recover Minkowski 
spacetime results by Wick rotating the final results. In Euclidean space, QCD's classical 
action in a covariant gauge dA a + a b a = takes the form 

S = /c&OCqcd + Asf), (1) 
where £qcd and £qf are given by 

£qcd = - F; v + igftT a -m q )% (2) 



jCqp = — < -r b a b a - b a (dA a ) + c a (dDc) a (3) 

and the notation is as follows. A ^ denotes the gauge field, ip q and ip q the quark fields, 
c a and c a the Faddeev-Popov ghosts, F" = d^A® — d v A a + gf abc A l \A C V the field strength 
and D^f = 5 ac + gf abc A b ^ the covariant derivative, g is the coupling constant, f abc are the 
structure constants of the gauge algebra, T a are the generators of the gauge algebra in the 
fundamental representation, m q are the fermion masses and a is the gauge-fixing parameter. 
The constants f abc and the generators T a are normalized so that f acd f bcd = N c 5 ab and 
tr(T a T b ) = \5 ab . SR regularizes QCD in two steps. First, it introduces a higher covariant 
derivative term and modifies the gauge-fixing part of the action so that now one has 

S A = J d 4 x (£ QCD + £ HC d + £' GF ) , (4) 



with 



£hcd = ^ (£>%,)" (D 2 Fn a (5) 



and 

£'gf = - f ^J^]^ ba - ba ^ + ^ ( dDc T ■ ( 6 ) 
Here A is a mass and f(d 2 /A 2 ) is a function that in momentum space is given by 
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Some simple power counting shows that the only superficially divergent 1PI Feynman di- 
agrams generated by S\ are the one-loop diagrams contributing to the two, three and 
four-point 1PI Green functions of the gauge field. Thus the modification of QCD's action 
along the lines of eqs. (f|)-(|6|) does not regularize the theory completely, but leaves some 
one-loop divergences unregularized. Before explaining Slavnov's idea to regularize the latter 
divergences, let us mention that the choice of f(d 2 /A 2 ) is somewhat arbitrary. Here we have 
chosen it as in (0) so as to ensure locality and make all a-dependent contributions finite by 
power counting. Taking e.g. f = 1 + (p 2 /A 2 ) also ensures locality but leaves a-dependent 
contributions unregularized at one loop (see ref. || for a discussion of this point). 

The second step in SR is to regularize the one- loop divergences generated by • Slavnov 
H proposed to do this by introducing Pauli-Villars determinants in the generating functional 
so that it reads 

z[j,x,c,c] = [ m m m m m $c] e -^ +s -^ 

Nf J, (8) 
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x n (det A t y ai/2 (det c t r n n (^^) Vqj 

i=l q=lj = l 



Here 

r \ Nl , - 

Source = / A J^A} + 6 V + E ( + 

is the source term coupling the fields A a , b a , ipq an d i>q to their external sources J aM , x a > C? 
and ( q , and the determinants det^lj, detCj and det^j are defined by 

(detA t )- 1/2 = J^A i \S{DA i ) ^v{~\ j ' J ' <&u AT{x)O l f v {x - y) A h f{y)) (9) 

det d = det (-L> 2 + M 2 ) (10) 
det T qi = det (ip + gA a T a + fM qj ) . (11) 

The parameters a, and rj q j are arbitrary real parameters satisfying the conditions 

I 

E^ + 1 = ( 12 ) 

Jq Jq 

E Vgj + 1 = m 2 q + E Vqj t4j = , (13) 

3=1 3=1 

Mi and \i q j are masses and the operator O i ^ / (x — y) in eq. (|9]) is given by 

OStix -y)= £A 7^%T7 r + M 2 5 ab g^5(x - y) . 
vj 8A a {x)8A h v {y) v ' 
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Strictly speaking, Slavnov only considered pure Yang-Mills theory, so he did not need to 
introduce determinants det T q j to regularize the divergences generated by quarks running 
along internal loops. We will see, anyway, that the determinants det T q j do not pose any 
problem and that the difficulties arise from the determinants det*4.j. It is very easy to see 
H that det^4j, detCj and detJF^ are gauge invariant. This, together with the gauge 
invariance of S\, implies that the functional Z[J, x, (, C] satisfies the same BRS identities 
as the unregularized functional one would construct starting from QCD's action S in eq. (|l|). 
It can also be shown || that conditions ( |T2"D and ([13|) ensure that Z[J, x, (,(] generates 
finite Green functions at one loop. So, all in all, one has a generating functional which is 
manifestly gauge invariant and that generates finite Green functions at one loop. 

It is very important to understand the regularization mechanism of one-loop divergences 
in Z[J, x, C>C]- To this end, let us consider the vacuum polarization tensor it generates. 
From the action S\, it receives the contributions of the diagrams in Fig. 1. In addi- 
tion there are the contributions from Ili(det Ai)~ ai ^ 2 (det Ci) ai Y\ q T7j (det Tqjf 113 in the mea- 
sure of the path integral. Using the very same techniques as for ordinary Pauli-Villars 
regularization of QED, it is very easy to see that conditions (|13|) imply that the product 
I\q I\j (det J-qj) Vq3 regularizes the divergences in diagram (la). As for diagrams (lb)-(ld), 
Slavnov []5] has argued using formal path integral manipulations that violate locality that 
the product ]li(det Ai)~ ai ^ 2 (det Ci) ai cancels the divergences in diagrams (lb)-(ld) provided 
eq. (|l2l) and the extra condition 

i=l 

are met. It has been shown || that, to check whether this is actually the case within 
the framework of local regularization, an extra regulator TZ must be introduced. It turns 
out that after introducing such a regulator and performing calculations at finite TZ, the 
2-point divergences that arise in ]li(det Ai)~ ai ^ 2 (det Ci) ai when TZ — ► cancel the di- 
vergences that arise in diagrams (lb)- (Id) when TZ — > 0, provided only condition ([TJD 
is satisfied. So indeed the sum of diagrams (la)-(ld) with the 2-point corrections from 
n,(det A)" ai/2 (detCi) Ql U q Uj (detJ r qj ) Vqj is finite if eqs. ([TJD and (0) hold, but to see it 



without giving up locality, an extra regulator TZ is needed. The fact that one has to intro- 
duce an extra regulator TZ means strictly speaking that by itself SR does not provide a local 
regularization of QCD. 

Having a generating functional that generates finite Green functions for finite values of 
the masses A, Mj and fi q j is not all. One has to devise a subtraction procedure that 
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removes the divergences associated to large values of the regulators A, Mj and /j, q j, while 
preserving gauge invariance. There have been claims in the past [[/J that such a procedure 
does not exist for pure Yang-Mills theory and that this is enough to kill SR. In ref. ||, 
however, it has been proved that this is not the case and the most general subtraction 
procedure consistent with gauge invariance for pure Yang-Mills theory has been given. Its 
generalization to QCD being straightforward, we will not present here. In what follows, we 
show that the functional Z[J, x,(i C] generates unphysical contributions that, after Wick 
rotation to Minkowski spacetime, spoil unitarity. 

3 Violation of unitarity. Let us now come to Minkowski spacetime, the correct framework 
to discuss unitarity. As is well known, unitarity implies that the transition amplitude Tfi 
for a physical process \i) — > \ f) must satisfy the relation 

2 lmT fl = £ (2tt) 4 5\p n - Pi ) T*f T m , (14) 

n 

where the sum is extended over all physical intermediate states |n) connecting |z) with |/) 
and p n denotes the momentum of the state \n). Consider the process fermion, antifermion 



going to fermion, antifermion. For this process, eq. (|1J) takes at first order in perturbation 
theory the form 

2ImT ljJf _ ff = £(27r) 4 5 4 (p n - Pi ) |^_ rf | 2 . (15) 

n 

The renormalized amplitude 7\ = T x ff^ff can be computed by first regularizing and then by 
subtracting the divergences associated to the particular regulator one has used. This way, the 
left-hand side becomes regularization and subtraction-dependent. Actually, regularization- 
dependent only, since different admissible subtraction schemes differ by finite local renor- 
malizations and these carry finite local radiative corrections which do not reach the imagi- 
nary part of the transition amplitude. The right-hand side, however, is regularization and 
subtraction-independent, since it only involves the Feynman rules of the unregularized the- 
ory. Hence eq. ([15]) can be viewed as a necessary condition that the particular regularization 
and subtraction prescriptions used to renormalize the theory must satisfy in order to pre- 
serve unitarity. The idea of our proof of violation of unitarity by SR is to compute Im T\ in 
any SR-based renormalization scheme and see that it does not satisfy eq. (|T3J). Now, since 
dimensional regularization (DR) preserves unitarity, the right-hand side in eq. (|i~5| ) is equal 
to twice the imaginary part of T\ as computed in any DR-based renormalization scheme. 
Thus eq. ( |T5D can be replaced with 

2ImT liSR = 2ImT liDR , (16) 
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where the subscripts SR and DR refer to the regularization method used to compute Ti. In 
the sequel we show that eq. ([16]) does not hold. 

As is well known, the imaginary part of the amplitude 7\ receives contributions from 
the diagrams depicted in Fig. 2, where all external legs are on shell and the shadowed 
bubble stands for the vacuum polarization tensor at one loop. To compute the contribution 
of these diagrams to IulT^sr, we proceed as follows. We first calculate the renormalized 
contribution of every diagram to the amplitude T 1SR in Euclidean space; once we have done 
this, we Wick rotate to Minkowski spacetime; finally, we take the imaginary part. As concerns 
the technical aspects of this computation, we note that the calculation of the renormalized 
contribution of any of the diagrams involved requires computing its limit A,Mi,n q j — > oo 
and subtracting the divergences associated to this limit. The evaluation of such limit is 
tedious but straightforward if one uses the techniques developed in refs. [f|] and [[J. For 
simplicity, and since transition amplitudes are independent of the gauge fixing parameter a, 
we will work in the Feynman gauge a — 1. 

We start by looking at diagrams (2a) and (2b). If we amputate the external legs, we 
are left in both instances with a 1PI diagram whose expression in Euclidean space in the 
Feynman gauge is 

N f 

G 4 (pi,P2,p'i; A) = g 4 T a T b T b T a £ I q (A) , 



9=1 



where 



/.(A) = I* ^ <>■ + J- m J ™' &±* 7 a> 3 -P2+k, A) D^k, A) 

J (27r) 4 [ (p 1 + k) 1 + m l q \ [ (p[ + k) 1 + m l q \ 



and 



A f g^v k^k v 



fc 4 + A 4 V k 2 k^ + A\ 
is the gluon propagator for the action 5*a in eq. (^) with a = 1. Using the Lebesgue 
dominated convergence theorem, it is straightforward to see that the A —>■ oo limit of / g (A) 
is well defined and equal to 

Hm j (A \ = a [ d4k 7^ (jfi±fjz m «) A + t- m i) ^ ( ]7 ) 

a-oo qy } J^7ry[(p 1 + ky + mi}[{ P , l + ky + mi}(p 1 -p 2 +kyk 2 ' 1 ; 

Hence no infinite renormalization is necessary. This is no surprise since the unregularized 
1PI four-fermion vertex G4 at one loop is finite by power counting and is given by the right 
hand-side in eq. ([H]). This implies that the renormalized four- vertex G4 at one loop as 
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computed with SR agrees with the renormalized four- vertex G4 as computed with any other 
acceptable regularization method, and in particular with DR. It follows then, after Wick 
rotating to Minkowski space and putting the external legs on shell, that the contribution of 
diagrams (2a) and (2b) to the imaginary part of T\ in any SR-based renormalization scheme 
is the same as in any DR-based renormalization scheme: 

TmT*g = JmT*gi. (18) 



Next we move on to diagrams (2c) and (2d). They both have the 1PI diagram in Fig. (3) 
as subdiagram. For the latter subdiagram, SR gives in Euclidean space and in the Feynman 
gauge 



d 4 k 7a ty 2 + V ~ m i) 7m (A + $ ~ m g) lv ^xu 



2N C t{ J (2tt) 4 [ (p 2 + k) 2 + m 2 ] [ (p! + k) 2 + m 2 ] 



D Xv {k,k) . 



Using the techniques in ref. [|J to compute the large- A limit, and the results in ref. [10 to 
rewrite the contributions that do not vanish in this limit, we obtain 



ig 3 T a N { 



a 1 r rl /A^\ 

lim r*(pi,p2,A) = ^ — T \ja - In ^7 

A^oo ' " ' ) 8n 2 2N C f\ I 2 Km 2 



+ Vo 



+ Ff*(m q ,p 1 ,p 2 )}, (19) 



where vq is a numerical constant and the finite part is given by 



Ffj, n (jfiq iPii P2 ) = / dx dyy\n 



m„ 



D(x,y) 



+ C dx [ l dy m2ql>l - 2m<? +P2 ~ 2yP ^ + Wl ~ y ^ lvb ^ 2 ~ 
Jo Jo D(x,y) 



(20) 



Here D(x,y) and p M stand for 



D(x, y) = m; 2 + (p 2 - Pi) 2 x (1 - x) + y (1 - y) p 2 



and 



[1 — x) Pi + XP2 ■ 



We see that r^(pi,p 2 ,A) diverges as A — ► 00. To remove the divergence, we perform the 
most general subtraction compatible with gauge invariance and obtain 



ig 3 T a Nf 

r;, re >i,p 2 ,a)= — 2 — £ 



q=l 



7m 



(Tfl 2 \ I TP? \ 1 1 

^J + ^)J+i^K,Pl,P2)}, (21) 
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a being the renormalization mass scale and v {m 2 q /<j 2 ) a finite function that does not depend 
on the momenta and that is only restricted by BRS invariance (in a minimal scheme, it would 
be zero). Let us now recall what DR gives. In DR, instead of eq. (|T^) one has 



lim , r^Oi^e, v) 



ig- 



8tt 2 2N, 



'-- q=l kZ 



1 , (V 2 

-+ In — 

e \mi 



+ F^(m q , Pl ,p 2 ) f . 



where v is the dimensional regularization mass scale and vq is a constant different from 
that in eq. (19). After renormalization, one obtains the same expression as in eq. (pi]), 
modulo finite local radiative corrections. This implies, after Wick rotating to Minkowski 
space and replacing the subdiagram in Fig. (3) with its renormalized expression, that the 
imaginary part of diagrams (2c) and (2d) is the same for SR-based renormalization schemes 
as for DR-based renormalization schemes: 



T rr2c,2d T rri2c,2d 

ImT l,SR = ImT l,DR 



(22) 



Proceeding analogously, it is easy to see that the same holds true for diagrams (2e) and (2f): 

ImT^l^lmT^. (23) 



Let us finally look at diagram (2g), and let us concentrate on the one-loop vacuum 
polarization tensor Il^(fc) hidden in it. The latter is made of two contributions. First there 
is the contribution (we will denote it with a prime) from diagram (la) and the determinants 
detJFgj that regularize it. Its expression can be computed following the very same steps 
as for the vacuum polarization tensor in Pauli-Villars regularized QED [pl| . After some 
calculations, we get 



N f r J, 



16tt 2 3 £^ lh lq3 V 



3=1 



<i:i 



7T n 



(k^ku - k 2 g^) , (24) 



where h(x) is the function 

h(x) =4x+ (1 - 2x) v 7 ! +4x In 



+ 4x + 1 



V Vl + 4x - 

and 7Tq is a numerical constant. Then there is the contribution from diagrams (lb)-(ld) and 
from the Pauli-Villars determinants det^lj and detCj. This contribution (we will denote it 
with a double prime) is the same as for pure Yang-Mills theory, and its A, Mj — > oo limit 
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has been computed for arbitrary a in ref. 
that in the Feynman gauge 



]. Borrowing the results from there, we have 



n";> (fc, a, Mi 

where 



9 2 N ( 



16vr 2 



c gab 



A 2 In 



A 2 



B 2 a i ln 



i=l 



M? 



Tin 



Ao-Bo 
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(kpK - V) > ( 25 ) 



(26) 



and 7Tq is another numerical constant. The actual values of A 2 and B 2 depend on the way 
the masses A and Mj are sent to infinity. For example, sending A to infinity while keeping 
Mi finite, and taking in the result Mj — > oo gives different A 2 and i?2 as proceeding the 
other way around. The difference A 2 — B 2 is, however, independent of the path followed to 
approach A, M, — > oo and is always given by eq. (p6|) . Summing the contributions (0) and 
(f25|), performing the most general subtraction compatible with gauge invariance, and Wick 
rotating to Minkowski space, we obtain for the renormalized vacuum polarization tensor 

,2 



n 



ab 

/iv, ren 



(k,a) 
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n SR (fc» S ab {k^ky - k'^g^u) 



16tt 2 



where 



„ „, . lliV c , f k 2 \ 2 f / m 2 \ /m 2 \ /m 2 

- v K-^) - 3 S IK- * ) + ln u) + <^ 



(27) 



(28) 



a is the subtraction point and Tr(m 2 /a 2 ) is an arbitrary real function carrying local finite 
radiative corrections constrained only by BRS invariance. For momenta k^ such that k 2 > 0, 
the vacuum polarization tensor picks an imaginary part since 



Imn SR (fc 2 ,a) 



lliVc 
6 



6{k 2 



n N{ , 2m 2 



9=1 



k 2 



/ i _ 4m 2 



A' 2 



(29) 



Let us compare this with the DR result. We recall that in any DR-based subtraction scheme, 
one has in the Feynman gauge that 



Im n DR (P, a) = ^ 9(k 2 ) + 9(k 2 - 4m 

6 6 g=l 



2m 2 



1 



Am 2 



(30) 



We see that the coefficient in front of 8(k 2 ) is different from that in eq. (E29|). As explained 
in ref. [f|, the difference ^ — | = | is originated by the Pauli-Villars determinants detAi- 
It is plain now that the imaginary part of the renormalized contribution 



2.y 



v(p 2 )rT a u{ Pl )] U(k 2 , a) 1 {k,k p - k 2 gj [u(p[)YT a v(p[) 
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of diagram (2g) to the amplitude T\ is not the same for SR-based renormalization schemes 
as for DR-based schemes. In other words, 

ImT^ImT^V. (31) 
Putting together eqs. QT|), fl2|), (|23|) and (0), we have that 



ImTi jSR ^ ImTi iDR , 

as we wanted to prove. 

4- Conclusion. At this point we can draw the following conclusions: 

(i) The regularization method proposed by Slavnov violates unitarity, the violation being 
produced by the Pauli-Villars determinants det Ai that the method chooses. Let us try to 
gain some intuition of why this is so. Assume that we naively switch off the regulators in the 
regularized path integral in eq. (|8|), that is to say, that we send the masses A, Mi and fi q j 
to infinity. Then we should recover the unregularized QCD path integral. However, this is 
not the case [g [|l2j. To see the latter, we rescale pf the Pauli-Villars field Af M" 1 A" 



in (det Ai) 1//2 , take the limit A, Mi — > oo, exponentiate S(DAi) and integrate over dfhc 
once OBby parts. This leaves us with 

(det Ai)- 1 ' 2 ~ j HAi] Hbi] exp | - - J d 4 x (A 2 + 2 ADbi) J as A, -»• oo . 

Completing the square in the exponent and performing the integral yields (det D 2 )- 1 ' 2 . 
Since each (det Ai) -1 ' 2 is exponentiated to the power ai and the a^'s satisfy eq. (fT2]), we 
obtain a factor (det D 2 ) 1 ' 2 . As for the determinants detCj and detJ r q j, it is straightforward 
to see that their naive limits Mi — > oo and fi q j — > oo give unity. Thus taking the naive 
A, Mj,/iqj — > oo limit in Z[J, x, (,(] yields the unregularized QCD path integral plus an 
extra (det D 2 ) 1 ' 2 . This extra determinant introduces propagating degrees of freedom that 
couple to the gluon field through the covariant derivative and which are not present in QCD's 
action. In other words, SR modifies QCD even at the tree level. Obviously the properties 
of the modified QCD are not the same as those of the true QCD. In the light of this, it is 
very easy to understand SR's violation of unitarity. What T^sr is really standing for is the 
transition amplitude T\ for the modified theory. By the cutting rules of 't Hooft and Veltman 
fl3|| , Im T 1) sr will receive contributions from the new propagating degrees of freedom. Hence 
there is no way IulT^sr will agree with the imaginary part of T^qcd = T^dr- All this 
discards SR as an acceptable regularization method. 
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(ii) Note that the diagrams whose regularization only involves A, namely diagrams (2a) to 
(2f), give the correct contribution to Ti qqd i n the A — ► oo limit. This indicates that the 
higher covariant derivatives terms in eq. (|5|) by themselves do not cause problems, in agree- 
ment with [[U]. The question that remains open is to supplement higher covariant derivatives 
with a suitable local regularization that preserves gauge invariance manifestly. Let us recall 
in this regard that for a local regularization method to be such, it must provide integrals over 
loop momenta which are finite by power counting (this is what local regularization is about). 
If a prescription does not provide this finiteness by power counting, it should not be called 
a local regularization; not even in the event that the various divergent contributions from 
different divergent Feynman integrals cancel among themselves when the latter are properly 
defined through yet another regularization. 
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Figure 1: Contributions from S\ to the vacuum polarization tensor. 
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